Bounding the dimension of bipartite quantum systems 
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Let us consider the set of joint quantum correlations arising from two-outcome local measurements 
on a bipartite quantum system. We prove that no finite dimension is sufficient to generate all 
these sets. We approach the problem in two different ways by constructing explicit examples for 
every dimension d, which demonstrates that there exist bipartite correlations that necessitate d- 
dimensional local quantum systems in order to generate them. We also show that at least 10 
two-outcome measurements must be carried out by the two parties altogether so as to generate 
bipartite joint correlations not achievable by two-dimensional local systems. The smallest explicit 
example we found involves 11 settings. 

PACS numbers: 03.65.Ud, 03.67.-a 



I. INTRODUCTION 

We consider the following standard bipartite measure- 
ment scenario Two distant, separated parties, con- 
ventionally called Alice and Bob, share a joint quantum 
state p. The two participants carry out measurements 
on their respective states and may obtain two possible 
classical outcomes ±1. Let us specify the situation, that 
Alice (Bob) has number of ±l-valued observ- 

ables Ai, . . . , ■ . ■ , Bmg)- In a particular run of 

the experiment Alice and Bob each measures one observ- 
able, Ai and Bj, getting the respective outputs G ±1 
and bj € ±1 and then makes the product Uibj. By re- 
peating this process many times they can form the joint 
correlation 



E{a,,bj) = {A.,Bj) = Jr{pA^ ® Bj 



(1) 



which is the expected value of the product atbj. These 
joint correlations for all {i,j} form an mArriB dimen- 
sional vector, designated by c. 

We are interested in the following problem. We are 
given a vector c of joint correlations E{ai,bj), with 
1 = 1,..., TUA and j = 1, . . . , rriB, for a particular number 
of settings ttia and m b ■ What is the minimum Hilbert- 
space dimension which is needed to reproduce this vector 
c? Is it the case that all sets of these vectors for any 
niA and tub are possibly reproduced by two-dimensional 
local Hilbert spaces? In a recent work wc provided coun- 
terexamples for small values of niA and itib (numerically 
for = 4 and m-A = 8; an exact treatment was carried 
out for TUB = 6 and tua = 15), which shows that indeed 
higher than two-dimensional local Hilbert spaces are re- 
quired in order to produce all the set of such vectors 3. 
The same conclusion, based on the gap between a lower 
bound on the Grothendieck constant in infinite dimen- 
sions [J and an upper bound on the three-dimensional 



* Electronic address: 
t Electronic address: 



t vertesi@dtp . at omki . hu] 



kfpal@atomki.hu| 



Grothendieck constant Q , has been obtained by Brunner 
et al. 0] ■ In this paper the idea of dimension witnesses 
has been put forward as well. Brunner et al. defined a 
d-dimensional witness as a certain kind of generalization 
of Bell inequalities [l|, such that all quantum correla- 
tions arising from observables on d-dimensional compo- 
nent Hilbert spaces satisfy the inequality. Hence, corre- 
lations which violate the inequality could be established 
only by measuring systems of dimension larger than d. 
Thus dimension witness is a useful tool for measuring 
the dimension of multipartite quantum systems relying 
only on experimental data. 

The results in i and i prove the existence of two- 
dimensional witnesses in bipartite systems based on only 
joint correlations. On the other hand, Perez-Garcia et 
al. gave a non-constructive proof of the existence of 
d-dimensional witnesses for tripartite systems for any di- 
mension d. In particular, they proved in Ref. Q in their 
Theorem 1 that for every dimension d, there exists a di- 
mension D, a pure state {ip) on (g) (g) C^, and a 
Bell inequality with ±l-valucd observables such that the 
violation by \^) is greater than \/d up to some universal 
constant. We also mention a recent paper by Wehner 
et al. 01, in which a method has been given to bound 
from below the dimension of a Hilbert space by relat- 
ing this problem to the construction of quantum random 
access codes. Their method works for quantum systems 
involving any number of parties (can be applied even for 
a single system) and is most useful for small number of 
settings. 

In this paper we tackle the problem of whether there 
exist dimension witnesses for bipartite systems built from 
two-outcome measurements for any dimension d. Note, 
that the tripartite case has been affirmatively answered 
in Q. On the other hand, the fact that this may also 
hold for the bipartite case has been conjectured in Ref. 
based on the plausible assumption that the Grothendieck 
constant 0] for infinite dimension is strictly larger than 
the Grothendieck constant for every finite order. In the 
numerical works of 0, we found various dimension 
witnesses up to dimensions of three and four, respec- 
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tively. using bipartite Bell inequalities up to five binary 
measurement settings = = 5. But there, be- 
sides joint correlations, local marginal terms were also 
involved. 

In particular, in Sec.[TT]of the present work, we discuss 
the link between joint correlations of two-party outcomes 
arising from two-outcome measurements and dot prod- 
ucts of unit vectors in the Euclidean space. Armed with 
this connection two families of Bell inequalities are pre- 
sented (in Sec, mil and in Sec. lIVp . proving the existence 
of bipartite correlations for any dimension d which re- 
quire at least d-dimensional systems in order to generate 
them. The first family (Sec. IIII[) involves a finite num- 
ber of settings for any finite d, while in case of the other 
family the number of settings is infinite. In Sec. Illli as 
a by-product, it is shown analytically that joint correla- 
tions of two parties resulting from two-outcome measure- 
ments, which require more than two-dimensional compo- 
nent spaces, cannot be generated by =4 and ~ 6 
or by less settings. That is, these numbers of settings are 
not enough to define a two-dimensional witness. On the 
contrary, examples are given for a two-dimensional wit- 
ness for settings itib — ^ and tua — 7 and for settings 
niB = 5 and uia = 6. The paper concludes with Sec. |Vl 

II. REPRESENTATION OF JOINT 
CORRELATIONS WITH DOT PRODUCTS 

This section contains two lemmas. The first one is bor- 
rowed from [ll| (their Lemma 2), establishing a link be- 
tween dot product of unit vectors in the Euclidean space 
and joint correlations originating from projective mea- 
surements in the Hilbert space, whereas the second one 
is a strengthening of this lemma regarding dimension wit- 
nesses. 

Here we state the first lemma without proof. 

Lemma 1. (Acin et al. flJl l) Suppose Alice and Bob 
measure observables A and B on a pure quantum state 
G (E) . Then we can associate a real unit vector 
a G M?"^ with A (independent of B), and a real unit 
vector b £ ffi^'^ with B (independent of A) such that 
E{a,b) ~ a-b. Moreover, if \^p) is maximally entangled, 
then we can assume that the vectors a and b lie in R ^ . 

Let us fix some notations regarding dimension wit- 
nesses [2f based on joint correlations. A d-dimensional 
witness [2|, in terms of joint correlations ([1]), is a linear 
function of these correlators, described by a vector m of 
real coeSicicnts Mij , such that 

rriA riiB 

m-c = ^^Af,j£;(a„6,) <I^'', (2) 

i=l 3 = 1 

for all correlations of form ^ with p in C'', and 
such that there are joint correlations for higher than d- 
dimcnsional systems for which m • c > W^. Note, that 



due to convexity arguments it is sufficient to consider 
only pure states p = in order to achieve the max- 

imum on the LHS of Q. 

Lemma [T] implies the following. Let us take an arbi- 
trary vector m in ([2]) defined by some coefficients Mij, 
and denote the maximum value 

Q'' = max{m- c} (3) 

c 

achievable by correlations ([1]) from (X) C^. Because 
of Lemma [T] we can associate with the elements of c, 
E{ai,bj), two normalized vectors di,bj e M" such that 
E{ai, bj) — di ■ bj. Defining 

niA ins 

S = ^5]Af,,a.-6„ (4) 

i=l 3 = 1 

and denoting the extremal value in the n-dimensional 
Euclidean space by 

T"= inax B. (5) 

Lemma ID implies that Q'^ < T" for n = 2d'^. 

This way we can bound the quantum value Q"^ achiev- 
able in d-dimensional component Hilbert spaces from 
above by dot product of unit vectors in M". On the 
other hand, due to Tsirelson's construction [l^, all dot 
products of unit vectors, a, 6 e M™, can be realized as 
±l-valued observables on a maximally entangled state 
on ®C°, where D = 2L"/2J. Then in light of this 
and the fact that Q'^ < T" for n = 2d^, constructing a 
Bell expression Q for which T" < T™, with n < m, im- 
plies a d-dimensional witness. The next Lemma sharpens 
this statement. 

Lemma 2. In the notation defined above by (0) and 
we have Q'^ < T" for n ^ 2d — \. Furthermore, if the 
Hilbert space is restricted to be real, then the above re- 
lation holds true with n = d. Thus, by constructing a 
Bell expression B for which T" < implies a d- 

dimensional witness with d ~ [^^^^J • 

In Sec. lIIII we present a family of Bell expressions, char- 
acterized by Bob's number of settings n = ms, for which 
one can prove, that the above condition T" < holds 
for the (n -I- l)th member of the family. In case of this 
family, however, we can give the gap (T"+^ — T") ex- 
plicitly only for small ji values. In Sec. IIVI another Bell 
expression is provided, where this gap can be calculated 
analytically for every n. 

Proof. Here we prove Lemma [2] and start by proving the 
relation Q'^ < T'^ in the case of Hilbert space ® 
restricted to be real. Since and B^ are ±l-valued ob- 
servables the normalization conditions {o.ki)'^ ^ ^ 
and X^zLi (^i;)^ — 1 ^^'^ satisfied for all k = l,...,d, 
where A^ = (a^i) and B^ = (5^^) are d x d symmet- 
ric matrices with real coefficients (designations i and 
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j are moved to upper indices for convenience). On 
the other hand, for given m, Q"^ = maxp {m • c} = 
max^^^. Mij{ijj\A^ ® B^lip), where maximization is over 

ah in C'' (g) and for aU observables A' and 

in C*. The state \tp) can be written in Schmidt form, 

IV-") = X^iLi "^il**); where tti can be chosen positive and 
their squares add up to unity. Then we have 

m • c ^ A'U,{^P\A' ^ B^i') = Yl ^kaiNuu (6) 

ij kl 

where Nu = E., + ak^lfc)/2- Note that in 

the case of real matrices and B^ the products a^^ib-^i 
and alf.bji^ are the same. However, this form will turn out 
to be useful in the proof of the complex case. Now let us 
flip the signs of a^; and ajj. for alH = 1, . . . , d in each pair 
{fc, /} if Nki is negative. By performing the necessary sign 
flips we obtain the matrices A^ defined by elements a).;. 
In this way the sum N^i ^ Mrj{hlily>^i + h\^b]^)/2 
becomes positive for any given {fc, V\ pair. Thus, as a 
result the value of m • c cannot decrease. Furthermore 
the normalization condition still holds for the sign-flipped 

matrices, Ef=i ('^m)^ 1- the following, through a 
chain of inequalities we obtain an upper bound for Q'^ . 
First of all, we can write 

^akaiNki < J2 {otl+aj)Nki/2 (7) 

kl kl 

= Y.alJ2 M^, {-alibi, + a?,6^,)/2 (8) 

kl ij 

= E"'E^«-^ (9) 

k ij 

< E E ^^^^^ • = E ^^'^'^ • 

k ij ij 

(10) 

where the unit vectors al,bl and a*,fr' are in the d- 
dimensional Euclidean space. In inequality ([7]) we used 
Nki < Nki and the relation 2xy < + holding for 
any real number x and y. Equality follows from the 
fact that the matrices A^ and B^ are symmetric. In equal- 
ity ^ we exploited that aj., and iP^i can be treated as the 
d real components of unit vectors = {a\i, . . . , 

and IP^ = ^1,2' ■ ■ ■ ' ^i-d)- -'^'^ ^^^t inequality we 
omitted index fc, but keeping in mind the normalization 
conditions for a|. and we have (? • a* = 1 and V -V = 1 . 
Finally the last equality owes to Efc "^fc^ = 1- According 
to ^ and ([5]), Q'' is the maximum of ctkCtiNki over 
all observables and states of dimension d. Thus, by use of 
the chain of inequalities (jlOp . we obtain the upper bound 

Q'^ < max V AUjo" ■ V = T'^, (11) 
where and bj are unit vectors. 



The proof regarding the complex case goes along the 
same line as in the real case, but in this case = [a].,) 

and B^ ~ {^*ki) ^'^^ definition Hermitian matrices 
and their components are complex valued, where * de- 
notes complex conjugation. Similarly as in the real case, 
d\i is defined by flipping the sign of a\i if Nki is nega- 
tive for a given pair {k,l). Then in the complex case, 
in place of the real number {alibi, + a\,.b\,.) one can 
write 2(Rea^;Re(6;^;) — Ima^.;Im(fo;^;)), and then in ^ we 
have dl = (Rea'fci, Im5'fci, . . . , d\,^, . . . , Redl^, Ima'fc^) and 
bi = {Rebl^,-Iuibl^, . . . , bi,^, . . . , Re6^^^, - W^J. Since 
components a^.^ and 6^^. are real, these vectors lie in the 
{2d — l)-dimensional Euclidean space, and they can be 
checked to be of unit length. Therefore, in the complex 
case in ([TT|) one has the unit vectors d^,V S R^**^^ and 
the RHS becomes T^**^^, which completes the proof. □ 

Some remarks about Lemma [2] are in order. First, 
provided that the real Hilbert-space result is given, one 
can also obtain the result concerning the complex case 
by mapping the d x d Hermitian observables to 2d x 2d 
real observables (as discussed in the Appendix of Q and 
in the multipartite setting in (l^). which latter matrices 
have the property that at each column there is at least 
one zero component, entailing the (2c? — l)-dimensional 
Euclidean vectors a* and b^ in the upper bound expres- 
sion on Q"^. Second, in light of the proofs in Refs. [3, Q 
regarding the qubit case, in Lemma [21 is equal to 
and in the respective cases of complex and real qubits. 



III. BOUNDS ON DIMENSIONS WITH FINITE 
NUMBER OF SETTINGS 

In this section an example for dimension witness is 
provided for every dimension d. Let us consider the 
expression B defined earlier in Sec. |TT]by (|4]). Accord- 
ing to Lemma [21 a Bell expression with a specific ma- 
trix m for which T" < (y denoting the maxi- 
mum value of (HI) in M") implies the existence of a d- 
dimensional Hilbert-space witness in a bipartite system 
with d= [^J. 

In Sec, nil A] it is shown that if niA < mB{mB — l)/2, 
then vectors from a smaller space than dimension ruB are 
enough to maximize B. Then in Sec. lIIIBl it is shown via 
a particular family of Bell inequalities (labeled by Bob's 
measurements 7tib) that when tua = 'rnB{mB — l)/2 4- 1 
for any given value of ms, the whole ms-dimensional 
Euclidean space is required to maximize B. This proves 
for every n = — 1 that T" < and implies that 

the Bell coefficients m = (My) of this particular fam- 
ily define a dimension witness in Hilbert-space dimen- 
sion d = L^^J 1 where tub is the number of settings on 
Bob's side. In Sec. IIII CI we determine the gap between 
the maximum value achievable in ttlb and in {rriB — 1)- 
dimensional Euclidean spaces for small tub, and we also 
discuss the question of minimal number of settings in or- 
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der for witnesses for two-dimensional component spaces 
to exist. 



A. Limits on tlie number of settings 

Let us start from expression (|4|). Without the loss of 
generality, we may suppose that mg < tua- The dimen- 
sions of the real vector spaces are large enough to accom- 
modate vectors that maximize the expression. From 



(12) 



it is clear that when B is maximal, each points to 
the same direction as the linear combination of the hj 
vectors it is multiplied with. Therefore, however large the 
number of Alice's measurement settings is, her vector 
space need not have a larger dimensionality than that of 
Bob. Bob's ms vectors can always be accommodated in 
an ms-dimensional space. What we will show here, is 
that if mA < ruBirnB — l)/2+ 1, then even vectors from 
a smaller space are enough to maximize B. 

Let B be the same as B with Alice's vectors a,; chosen 
optimally. B is only the function of Bob's vectors bj: 



Eh 

i=l 



(13) 



We can always find the maximum of B by maximizing B 
in terms of the unit vectors 6 of an m ^-dimensional vector 
space. The terms of the above equation can further be 
written as 



V ' i>3 



where 



Xij = bi ■ bj. 



(14) 



(15) 



The Xy, {j 



< i < tub) values determine the rela- 
tive directions of the unit vectors and, therefore, up 
to an irrelevant orthogonal transformation, the vectors 
themselves. We will regard ;S as a function of these 
N = iriBimB — l)/2 numbers. Besides |Xy | < 1, the 
variables must satisfy several other constraints. The do- 
main of the function is where the Gramian matrix, the 
symmetric matrix containing Xij and Xa = 1 is positive 
semidefinite. Its determinant, the Gramian is zero at the 
boundary of the domain, which means that the vectors bj 
are linearly dependent. If the Xij variables maximizing 
B lie there, then a space of less than niB dimensions is 
enough to accommodate the corresponding vectors. The 
whole ms-dimcnsional space is only required if the max- 
imum is not at the boundary of the domain. The yS as a 
function of Xij is non-analytical only where h = 0, but 



that occurs only at the boundary, as = means linear 
dependence of the vectors. 

From these it follows that if the vectors bj maximizing 
B span the ms-dimensional space, then all partial deriva- 
tives of the function must vanish there. Let us introduce 
a single index denoted by a Greek letter instead of the 
pair {i,j}, say let there be = {i — l)(i — 2)/2 + j. 
We also introduce the following notation: x,^ = Xij, 
Yku = MktMkj, and Ck = J2t ^^It- With this notation 



N 

2 ^ ^ ^ku^u 



At the maximum of B 

dB d 

dxn dx, 



dB d ""^ 1 
7^ = ^E''^ ""E^*^";" 



fe=i 



k=l 



(16) 



(17) 



must hold. This is a set of N homogenous linear equa- 
tions for 1/Zfc. If the number of variables niA is less than 
A'^-l- 1, the equations may only have nontrivial solution if 
no more than ruA — 1 of them are independent. If such 
a solution exists, it may define Ik only up to a constant 
factor. Therefore, Ik = y/^k, where Afe is a nontrivial so- 
lution. All Xk must be nonzero and must have the same 
sign, because Ik must be positive. If these conditions 
hold, from Eq. (|16p we get a set of ttia linear equations 
for and x^: 



N 



2 ^ Yki^x^ - y^/^k = -Ck 



(18) 



The solution of these equations can only correspond to a 
solution of our problem, if is positive, x^, are within 
the boundary of the domain and if it defines a maximum. 
However, if all these are true, the solution is not unique 
if the number of equations ttia is less than the number 
of variables + 1 = mB{mB — l)/2 + 1. In this case 
the equations are satisfied on a whole subspace of the 
space containing the domain of B. The function is con- 
stant in that subspace up to the boundary of the domain. 
Therefore, where the subspace crosses the boundary, the 
function will still have its maximum value. This way B 
can be maximized in a less than ms-dimensional space. 

To summarize the argument, if ttia is not large enough, 
the equations requiring the partial derivatives to be zero 
usually have no solution corresponding to a maximum. 
In this case the function will take its maximum value 
with its variables at the boundary of their domain. With 
the Bell factors having some very special values the 
equations may be solvable, but then they cannot be in- 
dependent. Therefore they cannot define the variables 
in a unique way. The function will have the same value 
on a subset of its domain, extending to the boundary. 
This means that either there is no solution or there is an 
infinity of them. Only vectors spanning a less than tub- 
dimensional space will maximize B in the first case, and 
there will exist such vectors maximizing B in the second 
case. 
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B. Family of Bell expressions 



The second derivatives of equation (|20p are 



When niA = 77ib(tob — l)/2+l = then the whole 

TOB-dimensional space may be required to maximize B. 
An example of such a case is the following: 

rriB 

B= a^J{b^-bJ)+JaoY,h, (19) 

where the corresponding Mij is defined through Eq. ([5]), 
parametrized by m^, and 7 > 0. Note that setting 7 = 0, 
we obtain the Z„ family introduced in By choosing 
Alice's vectors optimally, 



N 



N 



mB + 2j2x^ (20) 



j/=i 



Here we used the single Greek letter notation for the 
{i,j} pair we introduced earlier. The partial derivatives 
of equation (^0)1 are: 



dB 

dx„ 



1 



7 



V2 - 2xc 



"13 + 2 Y.u=l 



(21) 



which are zero at 



27^ 



ruB 



272 + mB{mB - 1) 



(22) 



Thus at the extremal point all nondiagonal entries of the 
Gramian will be equal, x < 1, and the diagonal ones are 
all 1 (the bi are normalized). Let us consider a determi- 
nant of size k X k with all diagonal elements having the 
value of p and all nondiagonal ones are equal to q. Let 
us subtract the (i — l)th column from the ith one in turn 
for i = k, k— 1, . . . , 3, 2, in this order. Then each column, 
except for the unchanged first one, will contain {p — q) in 
the diagonal, {q — p) just above it, and zero everywhere 
else. Then add the ith row to the (i — l)th in turn for 
i = k, k — 1, ... ,3,2, in this order. This way only the 
diagonal element will be nonzero in each column but the 
first one. with a value of {p — q), while the first element 
of the first column will be the sum of all original entries 
in that column, that is, p+ {k— l)q. Therefore, the value 
of the determinant, which is now the product of all diag- 
onal entries, will simply be [p+{k—l)q]{p — q)''^^. With 
p = 1 and q = X > —l/{mB — 1), which is true if 7 > 0, 
the determinant is positive for any k < ms- Therefore, 
the Gramian is positive definite, which means that we 
are not at the boundary of the domain. This is the only 
extremal point, and at that point 



B = T' 



■niB 



^/j^ + mB{mB - l)/2. (23) 



dB 



dxadxp 



( 

(2 - 2xo)-^l'^io.ii - 7 + 2 ^ X, 



-3/2 



(24) 

At the extremal point all nondiagonal elements a 7^ /3 of 
— 1 times the matrix are equal, = 7[?71b + mBirriB — 
l)a;]~'^/^ > 0, and also all diagonal elements a = /3 are 
equal, p' = (7' + (2 — 2x)~^/^ > g'. From the values of the 
determinants of matrices with such entries given above, 
it follows that minus one times the second derivative ma- 
trix is positive definite at the extremum; therefore the ex- 
tremum is a true maximum. There is no other extremal 
point in the domain of the function, so this one has to 
be the absolute maximum of B. This proves T" < 
for Bell expression (jl9p for n = tob — 1, entailing a di- 
mension witness d = [m s /2J for every m b according to 
Lemma [2J 

The 7 — yJmB /2 is an interesting special value, in 
which case x ~ Q; that is, all bi vectors arc orthogonal to 
each other. Geometrically, we may consider bi as edges, 
and bi — bj as face diagonals of a hypercube, while ^i 
is a vector pointing toward one of its space diagonals. If 
we decrease 7, that is, the weight of this space diagonal in 
the expression to be maximized, the geometrical object 
will flatten along this direction, and at 7 = it will 
collapse to become an (m^ — l)-dimensional object. The 
corresponding Bell inequality has one less measurement 
settings for Alice, and it belongs to family Z in Ref. 0]. 
When 7 7^ 0, the one extra term is enough to prevent the 
collapse, and the whole ms-dimensional space is required 
to accommodate the optimum object. 

We can determine the classical limit as the maximum 
value of 



B = 



i>j 



Zj\+"f 



i=l 



(25) 



with Zi = ±1 for alH < m^. The expression is permuta- 
tion invariant; therefore it only depends on the number 
of Zi ~ +1 values. Let it be k. Then X)j>j \ ~ = 
2k{mB — k), because \zi — Zj\ — 2 when Zi 7^ Zj, which 
occurs k{mB ~k) times; otherwise it is zero. At the same 
time I ^ 1^" ^ ("^s — ^)| = |™B — 2k\, therefore, 

B = 2k{mB — k) + 7|mB — 2k\. This expression has the 
same value at k and at tub — k, therefore it is enough 
to consider k < niB /2. It is easy to show that B takes 
its maximum at kmax, which is the non- negative integer 
nearest to {tub ~ 7)/2, and its maximum value, that is, 
the classical limit is Bd — {rn^ -I- 7^ — 4A2)/2, where 
A ~ \{mB — 7)/2 — kmax\- For 7 > thb this value is 
77713. 

As we have shown, to get the maximum value for this 
family of Bell inequalities we need rTis-dimensional vec- 
tor spaces. With (777,3 — l)-dimensional spaces we may 
only get a smaller value, but it is interesting to know 
how much smaller. In the following this gap will be de- 
termined for particular number of 7773 's. 
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C. Two-dimensional witnesses 

We have determined numerically the maxima achiev- 
able in (ms — l)-dimcnsional spaces for small ms cases. 
The solutions are nontrivial, they have different struc- 
tures for different ranges of 7. As three-dimensional 
spaces correspond to qubits 0, the Bell expres- 
sion ([TO)) defined above with tob = 4 and tua = 7 are 
especially interesting. They have the smallest number of 
measurement settings for one of the parties among corre- 
lation type Bell expressions (i.e., involving only joint cor- 
relations), whose maximum violation can not be obtained 
with qubits, while the other participant has just as few 
measurement settings as absolutely necessary. The max- 
imum value that can be achieved with three-dimensional 
vectors (i.e. with qubits when working with tensor prod- 
ucts of Hilbert spaces) as a function of 7 consists of three 
regions. For 7 = the four vectors bi point toward the 
vertices of a regular simplex. For finite 7 this simplex be- 
comes somewhat distorted, its shape will be a pyramid, 
whose base is a regular triangle and whose apex is above 
the center of the base. The maximum is given by such 
a solution up to 7 = 1.4153, a value very near to -\/2, 
where a less regular shape takes over. This is where the 
three-dimensional maximum value differs the most from 
the global, the four-dimensional one, their ratio here is 
1.0107161. In the third region, valid for larger 7, the four 
bi vectors point toward the corners of the square. 

Remembering (according to remarks at the end of 
Sec. ini) that real qubits correspond to two-dimensional 
Euclidean spaces and recalling the proof of Sec. IIII Al it 
follows that correlation Bell inequalities with m^, < 
3 do not require complex Hilbert spaces for their max- 
imal violation. On the other hand. Bell inequality (fTO)) 
with ms = 3 and tua = 4 and with 7=1 does require 
complex numbers to achieve maximum violation (in the 
Euclidean space Bob's optimal vectors form an orthogo- 
nal triad). Note, that this inequality is just the same as 
the 3x4 setting elegant Bell inequality introduced in [iJl ■ 
From Eq. it follows that the maximum value is 6 in 
M.^, corresponding to complex qubits. On the other hand, 
the maximum restricted to R^, i.e., to real qubits, can 
be obtained through semidcfinite programming . Ex- 
pression (|19p to be maximized can be brought to a form 
containing only quadratic variables subject to quadratic 
constraints (actually we have 14 variables with 7 con- 
straints). This non-convex problem can be solved via a 
series of convex relaxations of increasing size [l6| . This 
technique provides in each step a better upper bound to 
the global optimum. We obtain the upper bounds of 6 
in the first order and 5.8894 in the second order, which 
latter value coincides up to numerical precision with the 
value which can be attained numerically. We used the 
GloptiPol 3 package [l3| to solve this global optimiza- 
tion problem. This approach thus gives the gap between 
the maximum quantum value achievable in the complex 
and real qubit spaces through an exact treatment (the 
ratio being 6/5.8894 = 1.01878). The gap takes its max- 



imum at 7 = 1.3946 with a ratio of 1.0208047. On the 
other hand, using a heuristic method for the ms = 5 
and niA = 11, and for the = 6 and niA = 16 cases 
the maximum ratios are 1.0062317 (at 7 = 1.6396), and 
1.0041964 (at 7 = 1.7642), respectively. 

Now let us come back to Bell expression with 
ms = 4 and niA = 7. From our proof it does not follow 
that this is actually the one with the minimum total num- 
ber of measurement settings whose maximum violation 
requires higher-dimensional spaces than qubits. If we al- 
low both participants to have more than the minimum 
number of measurements, the sum may be decreased, 
with 1713= nT-A = 5 it is just 10. We have generated and 
checked numerically all 44685 nonequivalent nontrivial 
inequalities with My values restricted to 0, 1, and — 1, 
and we have found that qubits were enough to get the 
maximum value for each of them. We then allowed Afy 
to be 0, 1, —1, 2 and —2, while confining ourselves to sym- 
metric matrices. We found no case requiring more than 
qubits for maximum violation even among these 7.66 mil- 
lion inequivalent cases. Although this still does not prove 
that there is no niB = rnA = 5 correlation type Bell in- 
equality with this property, but it makes it likely. With 
ruA + 1^3 = 11 we found several examples even with al- 
lowing only 0; 1 and —1 values for My . With ms = 4 
and niA = 7 wc found 11 inequivalent cases (including 
the one discussed above with 7 = 1), and with ms = 5 
and rriA = 6 our extensive search gave 79 examples. Note 
that in these cases all of our results are due to numeri- 
cal search, since the complexity of this particular problem 
was too large to be handled by the semidefinite program- 
ming approach discussed previously. We are quite con- 
fident, though, about the results obtained by heuristic 
numerical computations. 

The necessity of four-dimensional Euclidean space in 
achieving the maximum value means that in the Hilbert 
space picture a pair of qubits is not enough 0, [l] . Ac- 
cording to the construct of Tsirelson [l^ (see paragraph 
above Lemma [3]), any value one may achieve in four- 
dimensional Euclidean spaces, one can also get with mea- 
surement operators in four-dimensional complex Hilbert 
spaces. We have calculated the maximum violation of 
the 79 TOs = 5 and tua = 6, the 11 = 4 and ttia = 7 
inequalities mentioned in the previous paragraph, the 7- 
dependent ms = 4 and itia = 7 case for a few 7-values, 
and for the = 4 and tua = 8 example intro- 
duced in Ref. Q with determining the appropriate mea- 
surement operators via numerical optimization according 
to Refs. Ud^. In all cases complex four-dimensional 
Hilbert spaces were required, smaller spaces were never 
enough. 



IV. BOUNDS ON DIMENSIONS WITH 
INFINITE NUMBER OF SETTINGS 

In this section we consider a bipartite Bell inequality 
with a continuum infinite number of two-outcome mea- 



7 



surement settings for each party, which will be proved 
to serve as dimension witnesses for arbitrary dimen- 
sions. Let the indices of the measurement settings be 
m-dimensional unit vectors, that is, elements of S"™"^, 
the surface of the unit sphere in i?™. Let the Bell co- 
efficients be proportional to the inner product of their 
indices, 



M(a;, y) = m{x, y) 



(26) 



where the factor of m has been introduced for the sake 
of convenience. By construction the coefficient matrix 
M{x,y) is positive semidefinite. Eq. (pS)) has been 
worked out by Grothcndicck [l^ in context of functional 
analysis. We use it in the following in matrix version for- 
malized in [l^ . The Bell expression with n-dimensional 
Euclidean vectors will then be: 



B: 



dcr{x) 



da{y)m{x, y)d{x) ■ biy), (27) 



where a denotes the normalized surface of the m- 
dimensional unit sphere S*™^^, a{x), {x £ 5"'^^) and 

^(y); (y G S"^~^) are the n-dimensional unit vectors cor- 
responding to the measurement operators of Alice and 
Bob, respectively. Like in the previous sections, we 
denote the scalar product of these vectors with a dot. 
Grothcndicck constructed this example to provide a lower 
bound for his constant Q. Actually, he calculated ex- 
pression ([27]) for a uniform distribution of x,y £ S""~^ 
for m oo but to our knowledge the optimality of the 
solution for any m and n has not been proved yet. Below 
we give an optimal solution for any value of m and n, 
and show that by fixing m oo the maximum value in 
R" is strictly increasing in n. As discussed in Lemma [51 
this will provide us a dimension witness for any finite 
dimension. The maximum value of the Bell expression 
achievable with n-dimensional vectors will then be: 



T" = sup m d(j{x) I d(T{y){x,y)b{y) 

(28) 

where d{x) has been chosen optimally, that is, parallel to 



h{x) 



da{y){x,y)b{y). 



(29) 



In the following an exact value will be given for the 
maximum value of this Bell expression for any dimen- 
sion n in the Euclidean space with respect to the clas- 
sical case a,b = ±1 (i.e., setting n = 1). Due to the 
lemma of Rietz (2^ this ratio can maximally be 7r/2 for 
a positive semidefinite matrix M. Due to Grothendieck 
[18| this ratio can be achieved for infinite m, n. This 
proves incidentally that the solution in infinite m, n is 
optimal. However, next we can calculate exactly the op- 
timal values T" for any n and m in contrast to the case 
in Sec. HITB] 

Due to the linearity of h{x), for each of its cornponents 
hi{x) there exist a, Zi £ S"^'-^ and a scalar [2l[ such 



that 

hi{x) = Aj(x, z,) (30) 
The hi{x) is maximal at z^, and its value there is A^: 

Ai = /li(zj) = / d(j{y){zi,y)bi{y). (31) 

From these it follows that the maximum value of B in R" 
may be written as: 



T" = sup m / da{x) 



Y.XUx,z,)^ (32) 



Let us introduce the following generalized spherical po- 
lar coordinates for the components of unit vector x: 

Xi = cos tpi 



Xi 



COS ipi smipfj_ {1 < i < m) (33) 

m — 1 

The integral of a function J-{x) on the normalized unit 
sphere is: 



f2lT 



/ d(T{x)J^{x) = / dip,n-i I dip,-, 
J s-^-^ Jo Jo 



difii sin™ ' ^ ipi ■ 



dipi sin™"^ ipiT 



2 Sm V3,„„2 . 

1 



2soSl . . . S,n~2 ' 

(34) 



where Si = sin* ifd^ = ^/^T{{i + l)/2)/r((i + 2)/2). 

Let us consider the case of n < m. Let us choose the 
basis such that the last n unit vectors e,„_„+i, em-n+2, 
. . . , e,„ span an n-dimensional subspace which contains 
all Zi. The value of the expression ([32]) does not change 
if we replace x in the integrand with its projection PnX 
onto this subspace. Let us define ||-Pna^|| ~ \J {PuX, Pnx) 
and introduce 



a unit vector in the n-dimensional subspace. Then 

da{x) 



(35) 



5]A|(a:',z,)2.||F„x|| 



T" = sup m 

(36) 

With our choice of the basis one can easily see from 
definition ((33)) of the generalized polar coordinates that 
||P„a;|| = sini^m_„ sin(p„i_„_i . . .sin(pi, depending only 
on the first m ~ n angles. At the same time, x' is inde- 
pendent of these variables, therefore, the two factors in 
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the integrand of equation (|36p may be integrated sepa- 
rately. The integral of the first factor will be an integral 
on S*"^^, while the value of the integral of the second 
factor, as one can easily verify by substituting T with 
||P„x|| in equation p4|) . is Sm_i/s„-_i. 

We will determine T" by constructing an upper limit, 
and then finding an explicit solution which saturates it. 
Let us consider 



Sk-l 



da{x) ^ 



(37) 



The function depends on b{y) (through and Zi), and 
T" = msm-i/s„-i supg(^)gg„_i Unn H n < 171, and T" = 
™™P&(y)GS"-i ^"i" otherwise. By using J f{x)dx < 



^ J dx J f'^{x)dx we get: 



Ukn < 



P n 1 " 



(38) 



To get the last expression the integral of {x, Zi)'^ on the 
/c-dimensional sphere has been carried out, which can 
easily be done using appropriate polar coordinates. For 
Er=i we get 



Up n 

^^^i^ j _ dcriy)'^Xi{zi,y)bi{y), 

A 1 J ^ ■ 1 



(39) 



by substituting one A^ factor in each term according to 
equation ((3T|) and by changing the order of the summa- 
tion and the integral. We may increase the integrand by 
replacing the vector b(y) by the unit vector parallel to 
h{y), that is the vector whose components are \i{zi,y). 
We may not necessarily be allowed to choose h{y) this 
way, as A^ and Zi are determined by b{y), and this choice 
may be inconsistent. However, we can overestimate the 
integral by making this replacement, therefore: 



< 



d(j{y)^ 



^A|(z,:,y)2. (40) 



Comparing the right hand side of this inequality with 
equation (|32p we can see that the maximum value it may 
take is nothing else but T" divided by m. Therefore, it 
follows that if n < m, 



T" < 



Sn-l 



n m 



that is, 



TO s: 



m — 1 



-^m— 1 



n St 



For the last form we used the identity nSnSn 



(41) 



(42) 



27r. 



Let us choose b{x) to be the normalized projection of 
X onto the n-dimensional subspace of the TO-dimensional 
space spanned by the last n members of the basis intro- 
duced earlier, that is 



6j(a;) 



~-m—n-\-ii • 



m — n-\-i— 1 

COS(pr,i-n+i Yl sil^'i^P- 

/J,— n+1 



(43) 



For the last form we used the form of generalized polar 
coordinates ([SS]) . and that ||Pna;|| = 11™=" sin V'a'- 
To get hi{ej) we have to integrate the product of 



cos iy9j Jl'^^^j^ sin iy9^ and bi{x) on 5™ ^ (see 



equation (HI])). From equation ([M]) one can see that for 
Lpj and for ^pm-n+i^ we have to integrate the product of 
the cosine function and a power of the sine function if 
j ^ m — n + i. Therefore hi{ej) = for these values of j. 
From this it follows that Zi ~ em-n+i, being hi{x) max- 
imal for that vector. Then we can explicitly calculate 
Ai from equation (|3ip by integrating {e„i-n+iT x)bi{x) = 

nm — n ■ -r-r^M — 1 . 2 ? nm~A 

p=i sm(p^n^=m-n-HiSin fycosi^i-n+t On S"" 

Substituting this into equation (|34p . performing 



the integrations and simplifying the expression by 

nr=m-„+i-|-l Sm-r-1, WC get 



A,: = 



nm — n 



7n — n-\-i— 1 



nm~-n-\-i 



^n — l Sn—i-\-l ^n—i 



^n — i — l 



nSn-1 



(44) 



For the last equality the identities Sn/sn~2 = {n — l)/n 
and 7is„s„_i = 27r have been used. Thus Xi is indepen- 
dent of i, as it has to be due to symmetries. Then by 
bringing A = A; in front of the integral in equation (|32p 



and using Zi 



-n+i we can sec that the integrand 



remaining is nothing else but ||P„a;||, whose integral is 
Sm-i/ Sn-i, as we have seen earlier. The result we get 
then for T" is equal to upper limit ([1^ . 

For n > TO, it is easy to verify that the upper limit is 1 
(instead of equation (|4T|) we get T" < a/T"), which is the 
well-known quantum limit for this inequality. It can be 
reached with TO-dimensional vector space, with b{x) = x. 
As the limit remains the same for n > m, there is no need 
for the extra dimensions. Although the Bell inequality 
involves an infinite number of measurement settings for 
each party, a finite dimensional space is enough to reach 
the quantum limit. 

Now we set to oo in (|42|) and calculate T" for dif- 
ferent n values. In this case T" = s„/s„_i. By choosing 
n = 1 we obtain the known classical limit 2/7r. Thus the 
maximum quantum violation (quantum limit per clas- 
sical limit) is (7r/2)s„/s„_i. For n = 2 the ratio is 
7r^/8 ^ 1.2337, corresponding to measurement on pairs 
of real qubits, while with n = 3, corresponding to com- 
plex qubits, it is 4/3. We have to go up to n = 5 to get 
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a maximum violation of 64/45 « 1.4222, larger than the 
value of \pi one can achieve with the CHSH inequality. 
For n — > cx) the maximum violation is the well-known 
7r/2, which is Stt/S « 1.1781 times larger than that wc 
can achieve with qubits. 

Most importantly the ratio s„/s„_i is a strictly in- 
creasing function of n. Notably for n even it is 

(V4) nri?~^ (2* + 1)V((2« + 1)^ - 1) and for n odd it 
is equal to (2/7r) Jll"!^^^^ (2«)V((2*)^ - !)■ Clearly both 
are strongly monotone functions of n. This entails T" < 
r£n+\ by ^ oQ for all 7j and owing to Lemma [5] proves 
the existence of dimension witnesses for any finite dimen- 
sion. 



V. CONCLUSION 

In the present work we focused our attention on joint 
quantum correlations arising from local measurements 
on bipartite systems. Tsirclson has established a con- 
nection between these joint correlations in Hilbert space 
and standard inner products of unit vectors in Euclidean 
space. Based on this result we give a proof that con- 
structing a correlation Bell expression for which T" < 
jin+i |-jjg vectorial picture (where T" denotes the max- 
imum value achievable in the n-dimensional Euclidean 
space) implies joint correlations which cannot be repro- 
duced in a d-dimensional Hilbert space, where d = [-^^^J . 
This defines a d-dimensional witness. For this sake, we 
discuss two particular families of Bell inequalities. The 
one in Sec. IIIII involves a finite number of measurement 
settings and = 771^(7715 — l)/24-l, whereas the 
other one in Sec. IIVI involves continuously many settings 
on both sides. Though in the former case we cannot give 
the difference j^'^+i — explicitly, except for small values 



of TTiB, in the latter case this difference can be analyt- 
ically calculated for all n. This conclusively proves the 
existence of dimension witnesses for arbitrary dimension 
d in a bipartite quantum system. Recently, we learned 
that Brief et al. arrived in Ref. [2^ at similar con- 
clusions. Besides, we discuss the minimum number of 
measurement settings arising in correlation Bell expres- 
sions in order to generate dimension witnesses. In this 
respect we prove that the numbers of settings = 4 
and niA = 6 (or less settings) are not sufficient to gener- 
ate a two-dimensional witness. In contrast, examples are 
given for a two-dimensional witness for settings = 4 
and TTiA = 7 and for settings tub = 5 and vtla = 6. It 
remains an open question whether for the pair of settings 
niB = 5 and niA = 5 there exists a two-dimensional di- 
mension witness or not. 

Due to Tsirclson's work [isj for a given finite num- 
ber of binary measurement settings, it is always possi- 
ble to generate all the joint bipartite quantum correla- 
tions of binary outcomes in a finite-dimensional Hilbert 
space. However, in generic Bell expressions involving 
local marginal terms as well, this may not be true. 
With respect to it, Navascues et al. [2^ asked recently 
whether there exist scenarios with finite number of mea- 
surement settings for both parties, for which all quantum 
correlations can be attained by measuring an infinite- 
dimensional entangled system. We leave this interesting 
problem as a challenge for future investigations. 
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